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A VARIATIONAL THEOREM FOR CREEP WITH APPLICATIONS TO PLATES AND COLUMNS 1 

By J. Lybll Sandebs, Jb., Habvey 6. McComb, Jb., and Floyd R. Schlechtb 


SUMMARY 

A variational theorem, is 'presented jor a body undergoing 
creep. Solutions to problems of the creep behavior of plates, 
columns, beams, and shells can be obtained by means of the 
direct methods of the calculus of variations in conjunction with 
the stated theorem. The application of the theorem is illus- 
trated for plates and columns by the solution of two sample 
problems. 

INTRODUCTION 

Interest in the various effects of creep in metals has 
intensified in recent years in the aeronautical field because 
of actual or envisaged operation of aircraft in an elevated- 
temperature environment. Much experimental work has 
been done to determine the creep stress-strain relations of 
various metals at elevated temperatures, and much data 
have been accumulated on the creep collapse of columns and 
plates. Analytically, creep problems for beams, col umn s, 
and plates are more difficult than the corresponding elastic 
problems, because the creep stress-strain relations (which 
generally include the possibility of elastic and plastic strains) 
are nonlinear. One consequence of this nonlinearity is that 
the distribution of the stresses through the thickness of the 
beam, plate, or column is not nearly so sim ple as in the 
elastic case. Mathematical analyses of these structural 
components have often been restricted to simplified models 
(H-section columns, for example) to avoid this difficulty. 

No known systematic procedures have been devised for 
reducing a given creep problem to a set of differential 
equations, although some efforts have been made to develop 
variational methods for use in attacking these problems 
mathematically. Hoff (refs. 1 and 2) has suggested the use 
of a minimum-complementary-energy principle based on an 
analog between the creep-strain rates and the elastic strains 
in a nonlinearly elastic body; however, elastic and plastic 
strains are neglected and the method does not apply readily 
to plate and column problems. Wang and Prager (ref. 3) 
have given two general extremum principles (for a body 
with elastic, plastic, thermal, and creep strains) analogous 
to the minimum-complementary-energy and minimum- 
potential-energy principles in elasticity. In the present 
report, an analog of Reissner’s variational theorem in 
elasticity (refs. 4 and 5) is formulated (for a body with 
elastic, plastic, and creep strains), and special forms of it 


are found to be convenient for handling creep problems of 
plates, columns, beams, and shells. The use of the vari- 
ational theorem presented herein is illustrated by applica- 
tions to two creep problems: creep collapse of columns and 
plate bending. 
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V 
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W 0 ,Wj 
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m 

n, 


V 
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r 


Young’s modulus 
integral to be varied (see eq. (7)) 
integrals defined by equations (46) 
second invariantof stress-deviator tensor 
radial and tangential bending moments, 
respectively, per unit length 
radial and tangential bending-moment 
coefficients defined by equations (44) 
compressive load 
buckling load for column 
surface of a body 

part of surface where stresses are pre- 
scribed 

part of surface where displacements are 
prescribed 
surface traction 

prescribed surface traction on S, 
displacement of neutral surface in ar-di- 
rection 

displacement coefficient in z-directiou 
(see eqs. (28)) 

displacement of neutral surface in y-di- 
rection; also used as volume integral 
deflection of neutral surface in 2 -direc- 
tion 

deflection coefficients (see eqs. (28) or 
(44)) 

radius of circular plate 
width of column 
thickness of column or plate 
length of column 
empirical constant in creep law 
unit external normal to surface of un- 
deformed body 

empirical constant in creep law 
lateral load intensity on circular plate 
radial coordinate 


1 Supersedes NAOA Technical Note 4003 by J. Lyell Sanders, Jr., Harvey G. McOomb, Jr., and Floyd B. Schlechte, 1957. 
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Sij 

stress-deviator tensor 

t 

time 

u 

displacement in z-direction 

u t 

displacement vector 

u { 

prescribed displacement vector on S d 

V 

displacement in y-direction 

w 

displacement in 2 -direction 

z,y,z 

coordinates 

7x«,7yi,7xf 

shear strains in zz-, yz-, and zy-direc- 
tions, respectively 

A 

amplitude of initial bow of column di- 
vided by column thickness 


Kronecker delta 

e 

strain, positive in tension 

«t J 

strain tensor 

t 

noncreep portion of strain tensor 

B 

*U 

creep strain tensor 

tr,tt 

radial and tangential strains, respec- 
tively 

e t> € t» e * 

strain in x-, y-, and 2 -directions, respec- 


e 

X 

M 

p=r/R 


tively 

angular coordinate 
empirical constant in creep law 
Poisson’s ratio 


cr stress, positive m tension 

<r {j stress tensor 

stresses defined in equation (27) 
vo,(Ti stress coefficients in equations (28) 

<T r ,<r t radial and tangential stresses, respec- 

tively 

Dots over quantities denote differentiation with respect 
to time, a single prime denotes the elastic-plastic part of the 
.strain, and a double prime denotes the creep part of the strain. 

FORMULATION OF GENERAL VARIATIONAL THEOREM 

For a body subject to creep, the equilibrium equations 
nnd strain-displacement relations are the same as those for 
an elastic body; the only difference in the mathematical 
formulation of a creep problem is in the stress-strain relation. 
In the present development, the nonlinear strain-displace- 
ment relations and equilibrium equations are used because 
of the intended application of the variational theorem to 
collapse problems of columns and plates. However, the 
strains are still assumed to be small, and the stress-strain re- 
lations are written as for infinitesimal strains. The total 
strain rate is separated into a creep part and a part which in- 
-cludes elastic and plastic effects. The creep part of the strain 
rate is allowed to depend on time and on the stress deviator 
but not on the stress rate or the first invariant of stress. In 
this report, the temperature of the body is assumed to be uni- 
form and constant in time. 

In tensor notation, a simple example of the creep law 
assumed is 

tij=f(Ji,t)s i} (1) 

where s {j is the stress deviator and Ji is the second invariant 
- StjSij. For purposes of the present report, the only restric- 
tion on the stress-strain relation for the elastic-plastic part 


of the strain is that the following relation hold (see ref. 6) : 

<’’ijde i j=ki j dcr t3 ( 2 ) 

The equation for the components of finite strain is 

(«<. j+Uj.t+Ut, (lit, t) (3) 

or, in differentiated form, 

«(j— 2 (“<. (a*, tUt' j) (4) 

The equilibrium equation is 

[(5«+'M’i.t)c r j*],j=0 (5) 

The equation giving the surface traction T, in terms of 
stresses and displacements on the boundary is 

( 6 ) 

where n } is the unit external normal to the surface of the 
undeformed body. 

The components of the tensor e i} are not true strains, and 
neither are the components of the tensor c r i} true stresses. 
However, the specific virtual work due to a virtual displace- 
ment biu is given by 

5W=ar ij 5e {j 

Therefore, the tensor <r {] may be considered to be the gener- 
alized force conjugate to the generalized displacement 
e u . In reference 7, this tensor is called k t) and in reference 
8 it is called o rf } . 

As in many of the variational theorems of the mathemati- 
cal theory of plasticity (see ref. 6), the variational theorem 
to be stated here is for stress rates and strain rates; that is, 
if the states of stress and strain throughout the body are 
known at a given instant, then application of the variational 
theorem singles out the stress rates or strain rates that 
actually occur from those rates of stress or strain considered 
in the enunciation of the theorem. Thus, the operator 5 
in the formal process of taking the variations is applied only 
to time derivatives of quantities and not to the quantities 
themselves. Because of this, for example, 56^=0 because 
depends upon the state of stress but not upon the stress 
rate. (See eq. (1).) 

The quantity to be varied was constructed by an inverse 
process guided by the form of Reissner’s theorem and tbe 
first theorem of Wang and Prager. By trial and error, the 
terms necessary to obtain the desired result were found for 

the integral to be varied. Let T t =T { be prescribed on the 
part S, of the bounding surface (S' of a body and let u t =u 
be prescribed on the remaining part S d of the boundary. 
The variational theorem to be proved is 51= 0 where 

^ — 2 (IV— 

f dS- f (u f —u t )T { dS (7) 

Js, J3 d 

In the first term of the volume integral, e {J is understood to 
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bo written in terms of displacements and velocities (eq. 
(4)). In the third term of the volume integral, is under- 
stood to be written in terms of stresses and stress rates, 
whereas i" {J is written in terms of stresses. 

The variation of 7 is given by 


57 = J" ^e t j8cr i} +& i j8e iJ +<r { /U twi 8ut, r 


1 

2 



dV- 


f TtSutdS- ( [{Ui-uJ&fi+fM,} dS ( 8 ) 
js, J s<i 

where 8e'j=0 as previously noted. The second term in the 
volume integral may be transformed by an integration by 
parts. The following formula applies: 

J ar { j8i t /lV= ^ 0’ijnj(8ik~\~'Uit,i)SUtdS — 

J* [ (S j) (T^] , j oiifd V (9) 

The third term may also be integrated by parts; thus, 


f <T i jU k ' i 8Uk, ] dV= f <r iJ n j u t , ( 8u>dS— f (cr ij u k ' i ), j 8u k dV ( 10 ) 

Jv Js Jv 


By using equations (9), ( 10 ), and ( 2 ), equation ( 8 ) becomes 


8 1 — J {tijScr tJ — [(Stj.+'M*; ,) 8 v.k — (<r, jfijt, t). j Silt . — 

(<U+««) W- i ]}dV + [<r {^(5(t+-Wt, t]8UtdS — 

( Tt&UidS— \ [(u t —u,)8'I\ J t-i\ 8u { ]dS ( 11 ) 

J S l J S d 


or 

8l~ J (e<j — 8u^^[{8 i t^Ut,{)<r{j]j‘^-dV + 

f (i' i —T i )8u t dS~ f (ui—u^sftdS ( 12 ) 

js, j s d 


The coefficient of each of the variations vanishes by virtue 
of the stress-strain relations, equilibrium equations, stress 
boundary conditions, and displacement boundary conditions, 
respectively. Thus, it is proved that 57=0 for the stress 
rates and strain rates that actually occur. 

A term may be added to the volume integral to include the 
effects of body forces if necessary. In some special cases, 
the boundary integrals should be modified if boundary condi- 
tions other than those on T t or u t are imposed. 

In the next section, a form for I appropriate for application 
to thin- plate and column problems is developed. 


THEORETICAL APPLICATION OF VARIATIONAL THEOREM 
TO PLATE AND COLUMN PROBLEMS 


The variational theorem proved in the preceding section 
provides a powerful means for deriving an appropriate two- 
D20507- — 00 0 


dimensional theory for thin plates subject to creep (or one- 
dimensional theory for beams and columns) analogous to the 
two-dimensional theory for thin elastic plates. If certain 
simplifying assumptions are made about the dependence of 
the various displacements and stresses on the 2 -coordinate 
(measured normal to the middle surface), the integrations 
with respect to 2 in the integral I may then be carried out. 
The equations resulting from 57=0 then involve only x, y, 
and t as independent variables. 

In this report, the same simplifying assumptions with 
regard to displacements are made that were made in the 
elastic case. In particular, the Bemoulli-Euler hypothesis 
is retained. This hypothesis can be regarded as a conse- 
quence of the equilibrium equations and the strain-displace- 
ment equations as applied to a thin plate. The transverse 
shear stresses vanish on the surface of the plate and are 
thus expected to be negligible (compared with inplane 
stresses) throughout the thickness. By any of the common 
stress-strain relations (elastic, plastic, or creep), it follows 
that the transverse shear strains are negligible. If w is 
assumed to be approximately independent of 2 , then from 


it follows that 


bin , bu _ 

^“s+sr - 0 

(13) 

bw , bv _ 

(14) 

%=U-z 

bx 

(15) 

v—V—z 5 — 
oy 

(16) 

8 

11 

(17) 


where U, V, and W can now be identified as the displace- 
ments of the middle surface. Equations (15), (16), and (17) 
express the Bernoulli-Euler hypothesis. 

A more rigorous treatment, in which the dependence of 
the various quantities on the 2 -coordinate is deduced by 
expanding them in powers of a thickness parameter, might 
be possible. In some derivations of plate equations, 77, V, 
and W are interpreted as weighted averages of u, v, and no. 
(See refs. 4 and 9, for example.) However, the consequences 
of assuming equations (15), (16), and (17) to be true are 
expected to lead to accurate results for thin plates. 

Equation (3) for the finite strains in terms of displacements 
is simplified by dropping all nonlinear terms except those 
containing powers or products of slopes. This simplification 
leads to the Yon K&rmfin equations for an elastic plate. 
The resulting equations for the strains are 


W_ dW 1 /bTYV 
dx 2 da* 2 V hr / 

(18) 

bV _bW , 1 /dTYY 

(19) 

by 2 by* 2 \ by ) 

W bV n mV dTYdTT 
*by + bx * jZ bxby bx by 

( 20 ) 
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e«=T„=TW=0 (21) 

In order to be consistent with, these approximations, only 
terms containing powers or products of rates of slope axe 
retained in the second term of the volume integral in I. 
(See eq. (7).) 

Before the integration with respect to z can be carried out 
in the expression for I, the dependence of the stresses upon z 
must be assumed. Unlike the elastic case, the stresses do 
not vary linearly through the thickness, even though the 
strains do vary linearly. The particular advantage of the 
generalized form of Reissner’s variational theorem for use in 
the present problem is that approximations for the stresses 
and strains can be made independently. Thus, there is no 
necessity to invert the stress-strain relations in order to 
determine appropriate approximations for the stresses; in 
fact, to invert the stress-strain relations would generally 
be impossible. 

Many possibilities exist for assuming the form of the 
stresses as functions of z, and each possibility leads to a 
different set of plate equations; thus, no particular system 
of equations can be set up which can be called the creep 
equations for a plate, column, or beam. Generally, it is 
advantageous to assume the form of the stresses and dis- 
placements as functions of x and y also, in which case all 
the integrations in I may be performed. Setting the variation 
of the resulting expression equal to zero leads to a set of 
ordinary differential equations in which the independent 
variable is the time. Further details of the application of the 
variational theorem to -creep-collapse and creep-bending 
problems of columns and plates can best be communicated 
by means of examples, two of which are given in the following 
sections. 

EXAMPLE OP APPLICATION OF VARIATIONAL THEOREM 
TO CREEP COLLAPSE OF COLUMNS 

A simplified treatment of the problem of creep collapse of 
the uniform, pin-ended, rectangular-section column shown 
in figure 1 is given as the first example. The plastic part of 
the strain is omitted, and the stress-strain relation is assumed 
to be 

( 22 ) 


where the second term on the right includes the effects of 
primary and secondary creep. The axial displacement u is 
assumed to bo given by 



Figure 1 . — Pin-ended rectangular-section column. 


and the axial strain e by 


bU 5 W , 1 /WY 

e dr Z dr 4 "^2 \ dr ) 


In the present case, the following terins enter into the 
expression for J: 




./W.bWb'W anU 


bx dr dr 2 

1 /dTFV 


dW \ 
dr s / 


2 atiaKiUt - i= 2 a \bi 


2 (*o+2eo)o'<i=2g+^ p ~U a m 

f T t UidS=—PU(l,t)=0 (P is constant) 

Js, 

•where Z7(0, f)= 0. The special form of I for the column is 
then 


-IKS 


. /dZ7 , bW bW d 2 


dr dr 


d *W \ 
dr 4 / 


(26) 


For simplicity, the variation of a in the z-direction is 
assumed to be linear; thus, 

< r =<r*+^ if (27) 

The variation of displacements and stresses in the r-direction 
is assumed as follows: 

W=hW 0 sin y ' 

U—U<ix 


<r °“E <r ° 

* P . TTX 
ai ~bh ai sm T 

J 

where h is the thickness of the column, b is the width, and 
l is the length. The dimensionless quantities on the right- 
hand sides of equations (28) are functions of time alone. 
The integrations in I may be carried out, and the result is 
(for m=3) 

I = WE (^oW-o+i ffili^o+StroWo 4 ^— I (jr 0 4 -f--^ oy)+ 
^ ciot/o (0 («*»*„+ 

| <W<r 0 +| cr/crijj (20) 


where P e is the buckling load of the column. The system of 
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differential equations to which the present approximate 
analysis leads is obtained by equating to zero the partial 
derivatives of I with respect to each of the dotted quantities; 
thus, the result is 


ScroTTo+g ( h+6 ,r ohto=0 

P ■ bhE ■ /P W 1 \ M 30 ) 

The solution of this set of equations depends on the initial 
conditions assigned. Suppose that the column has an initial 
bow in the form of a half sine wave before the load P is 
applied. If the load is assumed to be applied rapidly, but 
not rapidly enough to introduce inertia effects, then the 
initial conditions for the creep problem are those existing 
immediately after the load is applied. If the initial bow is 
given by 

, . . 7 TZ 

M sm — 

the initial conditions for the dimensionless quantities are 


TTo(0) = 


! — — 
Pc 


<r 0 (0) = -l 


The problem can now be induced to solving the following 
equation for W 0 : 

(j £ -l) W 0 -3Xpt^E(~J (w o +fw o ^=0 (32) 

Except for differences in notation, this equation is similar to 
one given by Kempner (eq. (27) of ref. 10) for an idealized 
two-element column problem solved by a collocation method, 
except that the constant % % is replaced by % and p is replaced 
by 1 . Equation (32) is easily integrated to give 


Pc 1 

p l 


1+fw 


l+fw 

W(0) 


The collapse time t 0 is found by taking the limit as W 0 -> » ; 
thus, 


^-1 
P 1 


jlog l+- 
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The foregoing analysis is intended merely as an example; 
more elaborate assumptions as to the stress variation through 
the thickness and the variation of the several unknowns 
along the length of the column lead to systems of nonlinekr 
differential equations which may be solved by some numer- 
ical method. 

EXAMPLE OF APPLICATION OF VARIATIONAL THEOREM TO CREEP 
BENDING OF PLATES 

Tor a second example, the bending of a simply supported 
circular plate under a constant uniform load is considered. 
(See fig. 2.) For simplicity, the analysis is restricted to the 
time during which the deflections are small so that the linear 
strain-displacement relations can apply. The general form 
for I in this case is 

I— 2 (35) 

Here, again, the effects of plasticity are neglected. The creep 
law is that given by equation (1) . The various terms in the 
integrand must be specialized to apply to the present plate 
problem. In polar coordinates the shear stress and shear 
strain are zero from symmetry, and the remaining stresses 



2 


Figure 2. — Simply supported circular plate with lateral load. 
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— erCr+ejO - 1 


(36) 


If a linear strain distribution through the thickness is 
assumed, 


. . d 2 !^ . z dli 7 . 

2 ^3T -^T a. 


(37) 


For the elastic part of the second term in equation (35) 

\ (ov4-o-< s — 2p<r r <rJ (38) 

For the creep part of the second term, 

(39) 

Further, it is assumed that the form of / is such that 

(40) 

In this example, the exponent m is taken to be unity; this 
assumption corresponds to the cubic law in the uniaxial stress 
case. Then, 


(‘IjO'ij—g X (<7>*+ o-, 2 — <r,<r,) [b r (2<r r — <r t ) + ar t (2cr, — <r T )] (41) 


The variation of the stresses through the thickness is 
assumed to be linear; thus, 


a t = 


12 M r z 
= h? ■ 

12 M,z 

= h 3 


(42) 


Doubtless, after creep is well established, the stress distribu- 
tion cannot actually be linear in the thickness direction, 
but, for purposes of this example, such an approximation is 
sufficient. When integration in the z-direction is carried 
out, 


/= - J7 £ { pMrW" +M t W' +~ p (M*+M t 2 — 

2 A M.) +~ ^(M r 2 +M t 2 -M r M,) [M r (2M- 


M,) +M, (2M, —M r )] 



(43) 


where B is the radius of the plate, p=r/B, and the primes 
denote differentiation with respect to p. The following 
expressions, based on the elastic solution for this problem, 
are assumed for the displacements and moments : 


W=h(l—p 2 ) (TFo+pWl) ' 

M T =jEh?M T 'o(l — p 2 ) ► 

M,=£%. ! (M, i0 +p 2 M,, 1 ) _ 


( 44 ) 


Equations (44) are substituted into equation (43) and the 
integrations carried out. Then, on equating to zoro tho 
partial derivatives with respect to the dotted quantities, tho 
following equations are obtained: 


■ftir, 0 _ f~2^f|, o +•3^1, 1=0 

3lbfr, i = 0 


W 0 +W 1 -2 (g) [2M,, 0 - m ( 3¥ w +* u )J 


=2 


144XSV7? 
5 V 


!) 7 . 


l^o-3 (|) 5 (2M,. 0 +M, 0 ) / 2 

3T^o+l^i — 6 (£) (3-M<.o+2-M (il — pM r .o) 


(46) 


where 


(M 2 ) 3 !^ JJpCI— p 2 ) (2 M r -M t ) (M*-M r M t +M,*)d P 

(Eh i yi 3 =£ P {2M t ~M T )(M/-M T M t +M*)dp 

(Eh*) 3 I 3 =£p 3 (2M t -M r )(M r i -M r M t +M t 1 )dp 


Y (46) 


Equations (45) can be solved for the five dotted quantities 
to give 

(|)* (37 ,_ / ,_ W ] 
*.=^(fj a 


3(l+p) 

8 


(3/3 I\ — 7j)J 


M r ,o=M,, 0 =-iM u =- g 


(3/3— /1—/3) 


y (47) 


where the expressions for I u I 3 , and J 3 are 


It — |q (3M|,o 3 3il/(,o J il7 (i 1 — !*■ — j\I t , j 3 ) 

I* = h (5M ‘- » 3 + 9M '. ° iM ‘. i+7M ti0 M ( . 1 2 +2i\f ( . , 3 ) 
J 3 =^ (14M,. 0 3 +32M,, 0 2 M t . 1 +27M (i0 M,. 1 2 +8M,, 1 3 ) 




^ (48) 


Equations (47) are a system of first-order nonlinear differ- 
ential equations with time as the independent variable. Tho 
initial conditions required to determine the solution are given 
by the elastic solution to the loaded plate. (See ref. 11.) 
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Thus, 


W( „,_KD‘3»o 

M,.(0)=M, a (0 

U 2 / -By l+3ju 

M.,( 0 )-- - 35 - 


(49) 


Since M r ,o=Mi,o and M ri o(0)=M tt o(fi), it follows that 
M r ,o=M ( , 0 . Equations (47) with initial conditions of equa- 
tions (49) can be integrated numerically by the modified 
Euler method (ref. 12). Numerical calculations were made 
for a plate having a ratio of radius to thickness of 50. The 
value of E used was 7.4 X 10 8 psi and the value of \E S used 
was 2.72 X 10 6 per hour. These values together with the 
cubic uniaxial creep law correspond approximately to char- 
acteristics of 2024-T3 aluminum alloy at a temperature of 
600° P. (See ref. 2.) Results showing the time history of 
the maximum deflection of a circular plate under lateral load 
are presented in figure 3. 


*Ns. 

/ 
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O 200 400 600 800 1,000 1,200 


Time, hr 

Future 3. — Time history of maximum deflection of simply supported 

R 

circular plate under lateral load. -^-=50; material, 2024-T3 alumi- 
num alloy at 600° F. 


CONCLUDING REMARKS 

A variational theorem for creep has been formulated which 
is an extension of a variational theorem developed by Reiss- 
ner. Various systems of equations leading to approximate 
solutions to problems of the creep behavior of plates, columns, 
beams, and shells may be obtained by using direct methods of 
the calculus of variations in conjunction with the stated 
theorem. The application of the theorem is illustrated for 
plates and columns by the solution of two sample problems. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 5, 1957. 
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